The Annals of Probability 

2007, Vol. 35, No. 5, 1998-2020 

DOI: 10.1214/009117906000001169 

(c) Institute of Mathematical Statistics, 2007 

S ■ STOCHASTIC DERIVATIVES FOR FRACTIONAL DIFFUSIONS 

o, 

S ' By Sebastien Darses and Ivan Nourdin 

Universite de Franche-Comte and Universite Paris 6 

^^ ' In this paper, we introduce some fundamental notions related to 

, the so-called stochastic derivatives with respect to a given u-field Q. 

OO ' In our framework, we recall well-known results about Markov-Wiener 

diffusions. We then focus mainly on the case where X is a fractional 
diffusion and where Q is the past, the future or the present of X. 

f^ ' We treat some crucial examples and our main result is the existence 



of stochastic derivatives with respect to the present of X when X 
solves a stochastic differential equation driven by a fractional Brow- 
i-Q ' nian motion with Hurst index H > 1/2. We give explicit formulas. 



Oh 



1. Introduction. There exist various ways to generalize the notion of dif- 
ferentiation on deterministic functions. We may think of fractional deriva- 
Tlj" ' tives or differentiation in the sense of the theory of distributions. In both 

^ , cases, we lose a dynamical or geometric interpretation of tangent vectors 

(velocities, e.g.). In the present work, we seek to construct derivatives on 
ff^ ' stochastic processes which conserve a dynamical meaning. Our goal is moti- 

■^ \ vated by the stochastic embedding of dynamical systems introduced in [2]. 

This procedure aims at comprehending the following question: how can we 



in 



o 






C^ • write an equation which contains the dynamical meaning of an initial or- 

dinary differential equation and which extends this dynamical meaning to 
stochastic processes? We refer to [3] for more details. 

Unfortunately, for most of the stochastic processes used in physical mod- 
els, the limit 

Zt+h — Zt 

C^ ' does not exist pathwise. What can be done to give a meaning to this limit? 

One of the main tools available is the "quantity of information" which we 
can use to calculate it, namely a given cr-field Q. The idea is that one can 
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2 S. DARSES AND I. NOURDIN 

remove the divergences which appear pathwise by averaging over a bundle 
of trajectories in the previous computation. This fact can be achieved by 
studying the behavior when h goes to zero of the conditional expectation 



E 



Zt+h — Zt 



h 



Q 



Such objects were introduced by Nelson in his dynamical theory of Brownian 
diffusion [9]. For a fixed time t, he calculates a forward (resp., backward) 
derivative with respect to a given a-field Vt which can be seen as the past 
of the process up to time t (resp., J^t, the future of the process after time t). 
The main class with which this can be done turns out to be that of Wiener 
diffusions. 

The purpose of this paper is, on one hand, to introduce notions which can 
be used to study the aforementioned quantities for general processes and, on 
the other hand, to treat some examples. We mainly study these notions on 
solutions of stochastic differential equations driven by a fractional Brownian 
motion with Hurst index H > 2- In particular, we recall results on Wiener 
diffusions (case H = ■^) in our framework. We prove that for a suitable a- 
algebra, the stochastic derivatives of a solution of the fractional stochastic 
differential equation exist and we are able to give explicit formulas. 

Our paper is organized as follows. In Section 2, we recall some now clas- 
sical facts on stochastic analysis for fractional Brownian motion. In Section 
3, we introduce the fundamental notions related to the so-called stochastic 
derivatives. In Section 4, we study stochastic derivatives of Nelson's type for 
fractional diffusions. We show in Section 5 that stochastic derivatives with 
respect to the present turn out to be adequate tools for treating fractional 
Brownian motion with H > 2- We also treat the more difficult case of a 
fractional diffusion. 

2. Basic notions for fractional Brownian motion. We briefly recall some 
basic facts concerning stochastic calculus with respect to a fractional Brown- 
ian motion; refer to [12] for further details. Let B = (5j)tG[o,T] be a fractional 
Brownian motion with Hurst parameter H € (0, 1) defined on a probability 
space {ft,J-',¥). This means that S is a centered Gaussian process with the 
covariance function 'E{BsBt) = i?//(s,t), where 

(1) RH{s,t) = Ut''' + s'''-\t "'2H 



s\ 



If if = 1/2, then B is a Brownian motion. From (1), one can easily see 
that E|i?t — i?<jp = \t — sp^, so B has a-Holder continuous paths for any 
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2.1. Space of deterministic integrands. We denote by £ the set of step 
M- valued functions on [0,T]. Let 7i be the Hilbert space defined as the 
closure of £■ with respect to the scalar product 

(l[0,t])l[0,s])H = RH{t,s). 

We denote by | • It^ the associated norm. The mapping l[o^i] i— > Bt can be ex- 
tended to an isometry between Ti and the Gaussian space Hi[B) associated 
with B. We denote this isometry by 99 i-^ B{ip). 

When /? G (2' 1)' ^^ follows from [14] that the elements of Ti may not be 
functions but distributions of negative order. It will be more convenient to 
work with a subspace of 7i which contains only functions. Such a space is 
the set 1 7^1 of all measurable functions / on [0,T] such that 

rp rp 

\f\f^^:=Hi2H-l)l^ 1^ \f{u)\\fiv)\\u-vr-'dudv<oc. 

We know that (1^1,1 ■ ||7^|) is a Banach space, but that i\'H\,{-,-)n) is not 
complete (see, e.g., [14]). 

Moreover, we have the inclusions 

(2) L\[0,T])cL^/^{[0,T])c\n\cn. 

2.2. Fractional operators. The covariance kernel RH{t,s) introduced in 
(1) can be written as 

rsAt 

RH{t,s)= I KHis,u)KH{t,u)du, 
Jo 

where KH{t,s) is the square integrable kernel defined by 

(3) KH{t,s) = CHs'/^~'' j\u-sf-^/^u''-^/^du, 0<s<t, 

where ch"^ = H{2H - l)/3(2 -2H,H - 1/2) ~^ and /? denotes the Beta func- 
tion. By convention, we set Knit, s) = if s > i. 

Let fC^ :£■ ^ L^([0,T]) be the linear operator defined by 

The following equality holds for any (pjip (^ £: 

{<P, ^)n = {JC*H^, IC*Hi^)LH[0,T]) = E(B((/.)i3(V')). 

JC^ then provides an isometry between the Hilbert space TC and a closed 
subspaceof L2([0,r]). 

The process W = {Wt)t£[o,T] defined by 

Wt = B{{JC*Hr\Mo,t])) 
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is a Wiener process and the process B has an integral representation of the 
form 

Bt= [ KH{t,s)dWs. 
Jo 

Hence, for any cp ^7i, 

B{cP) = W{IC*H<P). 

Let a, 6 G M, a <b. For any p>l, we denote by U* = LP{[a,b]) the usual 
Lebesgue space of functions on [a,b]. 

Let f € L^ and a > 0. The left-sided and right-sided fractional Riemann- 
Liouville integrals of / of order a are defined for almost all x a {a,b) by 

C/(x) = -i- /V - yr-'f{y) dy 



and 



lUi^) = y^ f\y - xT-'m dy, 



T{a) 

respectively, where F denotes the usual Euler function. These integrals ex- 
tend the classical integral of / when a = 1 . 

If / G la+iLP) [resp., / G 4°L(LP)] and a G (0, 1), then for almost all x G 
(a, 6), the left-sided and right-sided Riemann-Liouville derivative of / of 
order a are defined by 

and 

respectively, where a<x<b. 

We define the operator /Cjy on L^([0,T]) by 

{KHh){t)= f KHit,s)h{s)ds. 
Jo 

It is an isomorphism from L^([0,r]) onto /„+ (L^([0,T])) and it can be 
expressed as follows when H > ■^■. 

where h G L^([0, T]). The crucial point is that the functions of the space 
/q+ (L^([0,r])) are absolutely continuous when H > ^. For these func- 
tions (p, the inverse operator ICjj is given by 
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When iJ > |, we introduce the operator Oh on L^([0, T]) defined by 

(4) {OH^)is) := (|x;^) {^){s) = s^-'/^C's'/'-%{s). 

Let / : [a, 6] — > M be a-Holder continuous and g:[a,h\ ^ M be /^-Holder 
continuous with a + /3 > 1. Then for any s, t € [a, 6], the Young integral [18] 
!s f ^9 exists and we can express it in terms of fractional derivatives (see 
[19]): for any 7 € (1 — /3, a), we have 

(5) f fdg = {-ir f'D:^f{x)Dl:^gt.{x)dx, 

J s J s 

where gt-{x) = g{x) — g{t). In particular, we deduce that 



(6) Vs < t € [a, h] 



{f{r)-f{s))dg{T] 



< i^\f\o\g\f3\t 



\a+l3 



where k is a constant depending only on a,b,a and /?, and if h:[a,b] — > M 
and fi G (0, 1], then 

\h{t)-h{s)\ 
\h\^,= sup — -^. 

a<s<t<b \t-s\l^ 

2.3. Malliavin calculus. Let S be the set of all smooth cylindrical random 
variables, that is, which can be expressed as F = f{B{(j)i), . . . , B{(l)n)) where 
n > 1, / iM" ^ ]R is a smooth function with compact support and (pi (zTi. 
The Malliavin derivative of F with respect to B is the element of L'^{Q,'H) 
defined by 

'" df 
^sF = E JT^BicPi), . . .,B{(t>n))Us), s G [0,r]. 

In particular, DfBf = l[o,f](s). As usual, D-*^'^ denotes the closure of the set 
of smooth random variables with respect to the norm 

The Malliavin derivative D^ verifies the chain rule: if 93 : M" ^ M is C^ and 
(-Fi)i=iv," ^s ^ sequence of elements of B-^'^, then V3(Fi, . . . ,F„) € D^'^ and 
we have for any s € [0,T], 

Df (^(Fi, . . . , F„) = 5] -^(Fi, . . . , Fn)DfFi. 

The divergence operator 5^ is the adjoint of the derivative operator D^ . 
If a random variable u G L'^{'n,Ti.) belongs to the domain of the divergence 
operator, that is, if it verifies 

\E{D^F,u)n\ <c„||F||i2 for any F G 5, 
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then 5 {u) is defined by the duahty relationship 

E(F(5^(u))=E(D^F,u)7< 
for every FeB^'^. 

2.4. Pathwise integration with respect to B. li X = (^t)te[o,T] ^nd Z = 
{Zt)t£[Q,T] ai'e two continuous processes, then we define the forward integral 
of Z with respect to X, in the sense of Russo-Vallois, by 

(7) / ZsdXs = limucpe-i / Z,(X,+, - X,) ds, t G [0,r], 

provided the limit exists. Here "ucp" means "uniform convergence in prob- 
ability." If X (resp., Z) has a.s. Holder continuous paths of order a (resp., 
(3) with a + /? > 1 , then /q Zg dXg exists and coincides with the usual Young 
integral: see [15], Proposition 2.12. 

2.5. Stochastic differential equation driven by B. Here, we assume that 
H > 1/2. We denote by C^ the set of all functions whose derivatives from 
order 1 to order k are bounded. If a G C^ and b £ C^, then the equation 

(8) Xt = xo+ f a{Xs)dBs+ f b{Xs)ds, te[0,T], 

Jo Jo 

admits a unique solution X in the set of processes whose paths are Holder 
continuous of order a> 1 — H. Here, the integral with respect to B is in the 
sense of Russo-Vallois; see (7). Moreover, we have a Doss-Sussmann-type 
[5, 16] representation, 

Xt = <l)iAt,Bt), tG[o,r], 

where (f) and A are given, respectively, by 

- — {xi,X2)=a{(l){xi,X2)), (j){xi,0) = xi, xi,a;2GM 

OX2 

and 

A't = exp(-J ' a'{^{At,s))ds^b{(P{At,Bt)), Ao = xq, t e[0,T]. 
Using this representation, we can show that X belongs to 3^'"^ and that 
DfXt = a{Xs)expn b' {Xu)du + j^ a'{Xu)dB^ly^^t^{s), s,tG[0,r]; 
see [10], proof of Theorem B. 
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3. Notions related to stochastic derivatives. Let (-^t)fG[o,T] be a stochas- 
tic process defined on (0,J^, P). In the sequel, we always assume that for 
any tG[0,r], Zt e L'^{n,T,¥). For antG(0,r) and /i 7^ such that t + /i G 
(0,T), we set 

A ^ Zt+h - Zt 



h 
3.1. Stochastic derivatives in a strong sense. 

Definition 1. Set t G (0,r). We say that A^ (resp., B^) is a forward 
differentiating a- field (resp., backward differentiating a- field) for Z at t ii 
E[AhZt\A^] (resp., E[A_/iZ(|;B*]) converges in probability when hlO~^. In 
these cases, we define the so-called forward and backward derivatives, 

(9) DfZt=lmiE[AhZt\A% 

hiO+ 

(10) D^'Zt=limE[A^hZt\B']. 

hiO+ 

The set of all forward (resp., backward) differentiating cr-fields for Z at 
time t is denoted by M^ (resp., M^ ). The intuition we can have is 
that the larger A4^^^' is, the more regular Z is at time t. For instance, one 
obviously has that {0, Q.} € Mp^^ (resp., G A^^^*^) if and only if s 1-^ E{Zs) 
is right differentiable (resp., left differentiable) at time t. At the opposite 

extreme, one has that T S M^ (resp., G J^z ) ^^ ^'■^'^ only if s i—>- Zg is 
a.s. right differentiable (resp., left differentiable) at time t. 

Definition 2. We say that (^*,'B*)(g(o,T) is a differentiating collection 
of a- fields for Z if for any t G (0,T), A^ (resp., B*) is a forward (resp., 
backward) differentiating a-field for Z at t. If A^ = B* for any t G (0,T), we 
write, for simplicity, (>i*)tG(o,T) instead of (>i*,S*)jg(o,T)- 

On one hand, we may introduce the following definition. 

Definition 3. Set t G (0, T). We say that A^ (resp. B^) is a nondegen- 
erate forward a- field (resp., nondegenerate backward a- field) for Z at t if it 
is forward (resp., backward) differentiating at t and if 

(11) foranycGE F{DfZt = c)<l [resp., P{D^' Zt = c) <1]. 

For instance, if Z is a process such that s 1— > E(Zs) is differentiable at 
t G (0,T), then {0,0} is a forward and backward differentiating c-field at t 
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but is degenerate. Let us also note that condition (11) is obviously equiva- 
lent to \av{D^'Zt) / [resp., Var(D^*Zt) / 0] when D^' Zt G L'^i^) [resp., 
D^*Z4GL2(0)]. 

On the other hand, one could hope that such stochastic derivatives con- 
serve the property which holds for ordinary derivatives on functions: "it can 
discriminate the constants among the other processes." So we introduce the 
following. 

Definition 4. We say that (^*,S*)jg(o,T) is o. discriminating collection 
of a-fields for Z if (^*,;S*)ig(o,T) is a differentiating collection of u-fields for 
Z and if it satisfies the following property: 

(Vt e (0, T) , D^' Zt = D^'Zt = 0)^Z is a.s. a constant process on [0, T] . 

As in Definition 2, we write, for simplicity, {A^)te{o,T) instead of (^*, ;B*)fg(o,T) 
when A'' = S* for any t E (0, T). 

An obvious example of a discriminating collection of cr-fields for a process 
with differentiable paths is {A^ = J-^t ^ (0,T)}. If Z is a process such that 
s I— > E(Zs) is differentiable on (0,T), then the collection {A^ = {0,il},t € 
[0,r]} is differentiating, but, in general, not discriminating. 

Let us now consider a more advanced example. Let B = (-Bj)tg[o,T] be a 
fractional Brownian motion with Hurst index H G (1/2,1). Let us denote 
by Vt the a-field generated by Bg for < s < i and, if 5 : R — > M, by T^ the 
(T-field generated by g{Bt). 

Proposition 5. 1. For anyt G (0,T), Vt is not a forward differentiating 
a -field for B at t. 

2. For any even function 5 : M ^ M and for any t G (0, T), Tf is a forward 
and backward differentiating (but generate) a -field for Batt. 

3. For any t G (0, T), 7^ is a forward and backward differentiating and 
nondegenerate a -field for B at t. 

Proof. 1. We refer to Proposition 10 of [4] for a proof. This result is 
extended to the case of Volterra processes in the current paper; see Propo- 
sition 13. 

2. Since B and —B have the same law, it follows that 'E[lS.hBt\g{Bt)] = 
for any t G (0, T) and /i 7^ such that t + h^ (0) ^) • The conclusion follows 
easily. 

3. Using a linear Gaussian regression we can write 

n^HBtrn = (l + V^)^^-l-(VO^^ g^^ffgi -^ probability. 

2 ft— »o t 
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Since YaT{Ht~^Bt) > 0, 7^ is nondegenerate. D 

Thus, for fractional Brownian motion, stochastic derivatives with respect 
to the present (i.e., with respect to 7^"^) turn out to be adequate tools (see 
Section 5 below, for a more precise study). 

3.2. Stochastic derivatives in a weak sense. A way to weaken Definition 
1 is to consider stochastic derivatives as follows. 

Definition 6. Set t G (0, T) and let Ahea sub-a-field of JF. We say that 
Z is weak forward differentiable with respect to ^ at i if lim;jj^o+ E[yA/iZj] 
exists, for all random variables V belonging to a dense subspace of the closed 
subspace L^{n,A,F)cL'^{n,J^,F). 

We similarly define the notion of weak backward differentiation with re- 
spect to ^ at t by considering A^h^t instead of A^Zt- 

If the process Z is weak forward differentiable at t and such that the 
sequence {AiT_Zt)h is bounded in L^(fi), then we can associate with it a 
weak forward stochastic derivative with respect to A at t. Indeed, in that 
case, let us denote by Q the dense subspace involved. The linear form tp-.V >-^ 
lmii^^Q+ E[V A^Zt] defined on C L'^{Q,A,F) is continuous and so can be 
extended in a unique continuous linear form on L'^{^,A,¥), still denoted by 
i). Thus, there exists a unique Z^. G L^{n,A,F) such that tl^iy) = E[Z[V]. 
One can easily show that Z[ does not depend on 0. We will say that Z[ is 
the weak forward stochastic derivative of Z with respect to ^ at t. 

Remark 7. The boundedness of {AhZt)h in L'^{^) may appear to be 
quite a restrictive condition (e.g., it is not satisfied for a fractional Brownian 
motion). But it allows us to relate our notion to the usual notion of weak 
limit. 

If A^ (resp., S*) is a forward (resp., backward) differentiating cj-field for 
Z at t and if, moreover, the convergence (9) [resp., (10)] also holds in L^, 
then Z is weak forward (resp., backward) differentiable with respect to A^ 
(resp., B^) at t. But the converse is not true in general. 

Let T be the set of the so-called fractional diffusions X = (^t)tG[o,T] de- 
fined by 

(12) Xt = xo+ f asdBs+ f bsds, t£[0,T]. 

Jo Jo 

Here, a and b are processes which are adapted with respect to the natural 
filtration associated with B and X and satisfying the following conditions: 
o-gC" a.s. with a > 1 - i? and beL^{[0,T]) a.s. 
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Lemma 8. The decomposition ('i2) is unique. That is, if 

(13) xo+ / asdBs+ I hsds = XQ+ I asdBs+ / bgds, te[0,T], 
Jo Jo Jo Jo 

then xq = xq, a = a and b = b. 

Proof. The equality xq = xq is obvious and (13) is then equivalent to 



{as -ds)dBs 



which implies, by setting tk 

(K-^*Jl^t,+i-^tJ)^/^ 



kT 



{hs-hs)ds, tG[0,r], 



that 



•fe+i 



tk 



<C 



•fc+i 



hs)ds + 
., -b.)ds 



s — usju^a -r I yas - CFtk ) dBs + 

tk Jtk 

yn ftk+i 



•fc+1 



tk 



+ 



+ 



tk 

tk+1 



{as - at^ ) dBs 
{ds-atJdBs 



{as -a~t^)dB, 

l/H 
l/H. 



l/H 



We easily deduce, using (6), that 



n-l 



„l™o ^ 1^*^ - ^*^ l'^"" l^i'^+i - ^^>^\'^" = i^ probability. 

k=0 

But, on the other hand, it is easy to obtain (see, e.g., Theorem 4.4 in [7]) 
that 

n-l X 

J™oIIl^<fc-^*J^^'^l^*fc+i-^*fcl^^^= / Ws-as\^'"ds in probability. 
We deduce that a = a and so b = b. D 

In Section 4, we will see that the past of X € T before time t is not, in 
general, a forward differentiating cr-field at time t. We will see in Section 
5 that the present of X G T is, in general, a differentiating collection of 
(T-fields. 

However, X is weak differentiable for a large class of cr-fields. We intro- 
duce the set S'' of all r.v.'s (p{B{4>i), . . . , B{4>n)) G S such that (f>i,. .. ,(/)„_ are 
bounded functions. 

Let p be the set of all sub-cj-fields AcT such that L'^{n,A,F)nS^ is 
dense in L^(0,^, P). For instance, any a-field can be expressed as A'-^'^' = 
^{Bv,r <v<s) belongs to p (see, e.g., [12], page 24). 
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Proposition 9. Let Aep and t e (0,r). Let X eT be given by (VI), 
satisfying the following conditions: 

(i) the map si—>-bs is continuous from (0,T) into L^{Q,); 
(ii) for all s £ [0,T], as € B^'^ and sup^grgj^i E\Dfat\ < +00; 
(iii) E(|cr|g) < +00 for som,e p> 1 and a> 1 — H. 

Then X is weak forward and backward differentiable at t with respect to A. 

Proof. For simplicity, we only prove the forward case, the backward 
case being similar. Let t G (0,T). 
We write 

j-t+h j-t+h 

(14) Xt+h - Xt = at{Bt+h - Bt) + J bsds + J {as-at)dBs. 

First, we treat the second term of the right-hand side of (14). Let V G 
L^(r2, A, P) n S^. Since V is bounded and the map s 1— > 6s is continuous from 
(0,T) into L^{Qs), the function si— >E[l/6s] is continuous. We then deduce, 
by means of Fubini's theorem, that 



(15) 



lim— E 
h[0 h 



V 



t+h 



bsds 



ds = E[Vbt 



Then using inequality (6) and the hypothesis E(|cr|^) < +00, the following 
limit holds: 

1 r rt+h 
(16) ■• ■ 



lim — E 

h^O h 



ft+n 

Vj {as-at)dBs 



0. 



Finally, we show that the limit 



HO ^ * '' *^ 
exists. Since atV € D^'^ (see Exercise 1.2.13 in [11]), we have 
E[atiBt+h - Bt)V] = E[6''m^t+^)atV] 

= E[at{l[t,t+h] , D^'V)^] + E[Vm,t+^ , D^at)^]. 
Condition (ii) and the fact that V gS allow us, in particular, to write 
(17) E[at{Bt+h - Bt)V] = H{2H - l)(^t,h(^t, V) + ^t,,,(y, at)), 
where 

r pT pf.+h 

l^^^-^dyds 



^t^h{X,Y)=E 



T rt+h 

X I D?Y / \v 



/o Jt 

When X or y denotes at or V, Fubini's theorem yields 

rt+h 



^t,h{X,Y)=J^ f{v,X,Y) 



dv 
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with 

fiv,X,Y)= rE[XD^Y]\v-s\^''-''ds. 
Jo 

We have, due to condition (ii) and the fact that V ^ S , that 

\f{v,X,Y)-f{w,X,Y)\<C{X,Y) r\\v-s\''''-^-\w-s\''''-^\ds, 

Jo 

where C{X,Y) is a constant depending only on X and Y. 
The previous integral tends to when w tends to v. 
The continuity of the function v i— > f{v,X,Y) follows. 
Therefore, the limit 



exists and equals 
H{2H-1)E 



limh-'E[at{Bt+h-Bt)V] 
/i— >o 






U 
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4. Stochastic derivatives of Nelson's type. Let Z he a stochastic process 
defined on (i7,jr, P). We define the past of Z before time t, 

Pf :=?(Z,,0<s<t) 

and the future of Z after time t, 

Jf :=?(Z„t<s<r). 

If Vf and J-f are, respectively, forward and backward differentiating a- 

-pZ jrZ 

fields for Z at t, we call D^* Zt and DJ Zt, respectively, the forward and 
backward stochastic derivatives of Nelson^ s type in reference to Nelson's 
work [9]. In the sequel, we denote them by D^Zt and D_Zt, for simplicity. 

4.1. The case of Wiener diffusions. We denote by A the space of diffu- 
sion processes X satisfying the following conditions. 

1. X solves the stochastic differential equation 

(18) dXt = h{t,Xt)dt + a{t,Xt)dWt, Xo = xo, 

where xo eW^,b: [0, T] x M'^ ^ W^ and a : [0, T\xW^^W^® W^ are Borel 
measurable functions satisfying the following hypothesis: there exists a 
constant A' > such that for every x,y G M"^, we have 

sup(|(T(t,a;) -a{t,y)\ + \h{t,x) -b{t,y)\) <K\x-y\, 
t 

sup(|a(i,x)| + \b{t,x)\) < K{1 + l^l). 
t 
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2. For any t £ (0,T), Xt has a density pt- 

3. Setting Ojj = {aa*)ij, for any i,j € {l,...,n} any to ^ (Oi^) ^iid any 
bounded open set O C M'^, 

/ / |(}j(ajj(t,3;)pt(x))| dxdt < +CX3. 
Jto JO 

4. The functions b and (t,x) i— > — 4^ dj{aij{t,x)pt{x)) are bounded, belong to 

C^'^([0,T] X M'^) and have aU first- and second-order derivatives bounded 
[we use the usual convention that the term involving -^^ is ifpt{x) = 0]. 

These conditions are introduced in [8] and ensure the existence of a drift 
and a diffusion coefficient for the time-reversed process Xt ■= Xx-t- Follmer 
focuses in [6], Proposition 2.5, on the important relation between drifts and 
derivatives of Nelson's type. It allows the computation the drift of the time 
reversal of a Brownian diffusion with constant diffusion coefficient, both in 
the Markov and non-Markov case (see Theorem 3.10 and 4.7 in [6]). 

For a Markov diffusion with a rather general diffusion coefficient, we have 
the following result. 

Theorem 10. Let X e A be given by flSj. Then X is a Markov diffu- 
sion with respect to V^ and J-^ . Moreover, V^ and T-^ are differentiating 
and, in general, nondegenerate: 

DlXt = b{t,Xt), 
{D^Xt)i = bi{t,Xt) - ^—Y^dj{a,jit,Xt)pt{Xt)), 

where we use the convention that the term involving ]\ is if pt{x) =0. 

We refer to [3] for a proof; it is based on the proof of Proposition 4.1 in 
[17] and Theorem 2.3 in [8]. 

4.2. The case of fractional Brownian motion and Volterra processes. Let 
K be an L^-kernel, that is, a function iC:[0,T] x [0,r] — > M satisfying 
/[o T]2 K{t, s)"^ dtds < +00. We denote by ^^^ the right derivative of K with 
respect to t (with the convention that it equals +oo if it does not exist). 

We assume, moreover, that K is Volterra, that is, that it vanishes on 
{(t, s) € [0,T]^:s > t} and is nondegenerate. In other words, the family 
{K{t, •), t E [0,T]} is free and spans a vector space dense in L^([0,T]). With 
such a kernel K, we associate the so-called Volterra process 

(19) Gt= f K{t,s)dWs, 0<t<T, 

Jo 
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where W denotes a standard Brownian motion. The assumptions made on 
K imply, in particular, that the natural filtrations associated with W and 
G are the same (see, e.g., [1], Remark 3). 

Proposition 11. Let t e (0,T) and G be a Volterra process associated 
with a nondegenerate Volterra kernel K satisfying the condition 

The forward Nelson derivative D^Gt at t exists if and only if Jq q^ {t, s)^ ds < 
+0O. In this case, we have D^Gt = /q q^ {t,s)dWs and Vp is nondegen- 
erate at t if and only if Jq » (t, s)^ ds > 0. 

Proof. We adapt the proof of [4], Proposition 10. Using the represen- 
tation (19), we deduce that 

E[A^,Gt\V^]=E[AhGt\V^] 



I [ [K{t + h,s)- Kit, s)] dWs =: Zh. 
n Jo 



h 

Note that Z = {Zh)h>o is a centered Gaussian process. First, assume that 
/q ^-Q^{t,sy ds = +00. It is a classical result that if Z^ converges in proba- 
bility as ^ J, 0, then Var(Z/i) converges as h [0. But, from Fatou's lemma, 
we deduce that 

CX). 



limmfYav(Zh) > ^—(t,s)'^ds 
hio Jo at 



Thus, Zfi. does not converge in probability as /i | 0. Conversely, assume that 
/q g^ (t, s)'^ ds < -l-cx). In this case, assumption (20) implies that Z^ —^ 
Jq ^-Q^{t,s)dWs in probability as /i | 0. In other words, D^Gt exists and 
equals Jq ^^-Q^{t,s)dWs. We easily deduce that Vf is nondegenerate at t if 
and only if Jq ^(i, sf ds>0. D 

The result of Proposition 10 in [4] is then a particular case: if B denotes 
a fractional Brownian motion with Hurst index H G (0, 1/2) U (1/2, 1) and if 
t G (0, T), then D^Bt does not exist. Indeed, we have Bt = Jq Knit, s) dWg, 
where Kjj is the nondegenerate Volterra kernel given by (3) and verifying 



dt '' ' "\s 
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Remark 12. For a stochastic process Z, let us define 

(21) ^{Z)=Leb{t € [0,r], D^Zt exists}. 

For instance, if i? is a fractional Brownian motion with Hurst index H G 
(0, 1), then C{B) =T if H = 1/2 and ^{B) = otherwise. A real c € [0,T] 
being fixed, it is, in fact, not difficult, using Proposition 11, to construct a 
continuous process Z such that £,{Z) = c. For instance, we can consider the 
Volterra process associated with the Volterra kernel 



K{t,s) 



it- 
0, 



\Hit) 



with H{t) 



if s < t, 
otherwise, 

0, 

(t-c] 



A 1/4, 



if i < c, 
if t > c. 



The study of backward derivatives seems to be more difficult. Among the 
difficulties, we note that it is not easy to obtain backward representation of 
fractional diffusions (see [4]). However, for a fractional Brownian motion, we 
are able to prove the following proposition. 



Proposition 13. Set H > 1/2. The limit 

Bt — Bt-h 



limE 



h 



T, 



B 



exists neither as an element of U'{Q) for any p€ [l,C)o) nor as an almost 
sure limit. 



Proof. Fix t G (0,r) and set 
Bt - Bt-h 



G/.:=E 



h 



T 



and Zh := E 



Bt-B, 



t~h 



Bt,Bt+h 



Since {Gh)h>o is a family of Gaussian random variables, it suffices to prove 
that Var(G/j) diverges when h goes to 0. 

We have Z^ = E[Gh\Bt,Bt+h]- So, by Jensen's inequality, Z'^ < E[Gl\Bt, 
Bt+h] and Var(Z/j) < Var(G'/i). Let us show that lim/i^o Var(Z/j) = +oo. 

The covariance matrix of the Gaussian vector (Bt^h — Bt, Bt, Bt+h) is 



V* MJ ' 

where a = Vav{Bt-h -Bt),v = {R{t -h,t)- R{t, t) ■,R{t-h,t + h) - R{t,t + 
h)) and 

/ R{t,t) R{t,t + h) 



M: 



\R{t,t + h) R{t + h,t + h) 
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Since dh := i?(i, t)R{t + h,t + h)-R{t + h, tf / 0, M is invertible. There- 
fore, hZh = vM~^Q* , where Q = {Bt, Bt+h)- Since M = ^[Q*Q], we deduce 
that 

\s.T{hZh) =^[vM~^Q*{vM~^Q*)*]=vM-^v* . 

Hence, 

Yax{hZh) = 3-(i?(i + h,t + h)vl - 2R{t + h, t)viV2 + R{t, t)vl). 
dh 

This expression is homogeneous in t^ , so we henceforth work with t = 1. 
Tedious computation gives d^ ~ fi^^ as h [0. Moreover, we note that V2 = 
vi + c/i^^, where c is a constant depending only on H. Thus, 

dh VarihZh) = t;ic/i2^(l - (1 + hf^ + h^^) + h^^vf + c^/i^^ 
Since 2H > 1 and the function x i-^ x is derivable, the quantities ^ and 
i-(i+fe)2g^fe2g ^^^^^^g^ as /i i 0. But 2H <2 and ^^|^ = h^^'"^ -^ +oo as 
/i i 0. Thus, 

limVar(Z/j) = +cx), 

which concludes the proof. D 

4.3. The case of fractional diffusions. 

Proposition 14. Let X eT be given by (^12^ and satisfy the following 
conditions: E(/q \bs\ ds) < +oo and E{\a\^) < +oo for some p> 1 and a > 
1 — H . If for any t E (0,T), cr^^O a.s., then for almost all t S (0,T), V^ is 
not a forward differentiating a -field for X at t. 

Proof. Remember we assumed that a and b are adapted with respect 
to the natural filtration associated with B and X, see (12). In particular, 
we deduce from (12) that V^ C V^ . Since we can also write 

■ ds, 

Jo as Jo CTs 

we finally have Vi 

Thus, we deduce that ^[/\hBt\V^] =E[A?,Si|Pf] does not converge in 
probability as /i | 0, as a consequence of Proposition 10 in [4] or Proposi- 
tion 13 of the current paper. 

Consider expression (14). The hypothesis E f^ \bs\ds < +oo allows us 
to use the techniques of the proof of Proposition 2.5 in [6] to show that 
j^E[J^ bsdslVf] converges in probability for almost all t. Using inequality 
(6) and the hypothesis E(|(t|^) < +oo, we can finally conclude that V^ is 
not a forward differentiating cr-field for X at almost all times t. D 



Bt = 


/ -dXs- 
Jo OS 


Jo <y, 


Vf. 
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4.4. The case of fractional differential equations with analytic volatility. 

Proposition 15. Let X €E be given by (8) and let t G (0,r). We 
assume, moreover, that a is a real analytic function. Then V^ is a for- 
ward differentiating a-field for X at t if and only if a = 0. In this case, 
V = {"Pj , t G (0,r)} is a discriminating collection of a -fields and Vf is 
degenerate at any t G (0, T) . 

Proof. If o" = 0, then X is deterministic and differentiable in t. Conse- 
quently, Vf^ is a forward differentiating a-field, but is degenerate. Assume, 
now, that ct ^ 0. According to the Bouleau-Hirsch optimal criterion for frac- 
tional differential equations (see [10], Theorem B), we have that the law of 
Xt is absolutely continuous with respect to the Lebesgue measure for any 
t [indeed, we have inta~i({0}) = 0]. We deduce that P(<T(Xt) = 0) = for 
any t, since Leb(o"~^({0})) = (cr has only isolated zeros). Proposition 14 
allows us to conclude that V^ is not a forward differentiating a-field. D 

Remark 16. The case where a is not assumed to be analytic seems 
more difficult to handle. We conjecture, however, that in this case, V^^ is 
a forward differentiating u-field for X if and only ii t Kt^, where t^ is the 
deterministic time defined by 

t^ = inf{t > : Xt ^ int o-~^({0})} 

where {xt)t£[o,T] is the solution to xt = xq + /q b{xs) ds. If this conjecture is 
true, we would have S,{X) = t^; see (21). 

5. Stochastic derivatives with respect to the present. 

5.1. Definition. A consequence of Proposition 11 is that the a-field V^ 
generated by Xs,0 < s <t (the past oi X), cannot be used for differentiating 
when we work with fractional Brownian motion. Moreover, we stress the 
following important fact: the Markov property of a Wiener diffusion X £ A^ 
implies that taking expectations with respect to V^ produces the same effect 
as taking expectations only with respect to Xf. The following definition is 
then natural. 

Definition 17. Let Z = {Zt)te[o,T] be a stochastic process defined on 
a complete probability space (J7, JF, P) and for any t G (0,r), let 7^^ be the 
a-field generated by Zf. We say that Z admits a forward (resp., backward) 
stochastic derivative with respect to the present t G (0, T) if 7^^ is a for- 
ward (resp., backward) differentiating a-field for Z at t. In this case, we set 

DlZt := ofZt (resp., D^Zt := D^' Zt). 
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Example 18. Let i? be a fractional Brownian motion with Hurst index 
He (0,1) and tG(0,r). Then 

r Ht-^Bt, iiH> 1/2, 

DlBt =\q, iiH= 1/2, 

I, does not exist, ii H < 1/2, 

and 

r Ht~^Bt, iiH> 1/2, 

D^Zt = \t-^Bu \iH = 1/2, 

y does not exist, if i7 < 1/2, 

(see also Proposition 5). In particular, we would say that the fractional 
Brownian motion with Hurst index H > 1/2 is more regular than Brownian 
motion {H = 1/2) because of the equality between the forward and backward 
derivatives in the case H > 1/2, contrary to the case H = 1/2. We can 
identify the cause of these different regularities: the covariance function Rh 
is differentiable along the diagonal (t,t) in the case H > 1/2, while it is not 
when H=l/2. 

5.2. Case of fractional differential equations. We denote by H the set of 
fractional differential equations, that is, the subset of T whose elements are 
processes X = (^t)t6[o,T] solving (8) with cr £C^ and b gCI- 

In the sequel, we compute D^Xt for X G H and t G (0, T). Let us begin 
with a simple case. 

Proposition 19. Let X eE be given by (8) and let t G (0,r). Assume, 
moreover, that a and b are proportional. Then X admits a forward and a 
backward stochastic derivative with respect to the present t, given by 

(22) DlXt = DlXt = Ht~^a{Xt)Bt + b{Xt). 

In particular, the present T/^ is nondegenerate at t if and only if a{xQ) ^ 
and the collection of a -fields T^ = {T^^ , t G (0,T)} is discriminating for X . 

Proof. We will only provide the proof for D'^Xt, the computation for 
D_Xt being similar. Assume that b{x) = ra{x) with r G M. Then Xt = 
f{Bt + rt), where / : M ^ M is defined by /(O) = xq and /' = a{f). If a{xo) = 
0, then Xt = xq and D^Xt = = a{Xt)Ht-^Bt + b{Xt). If a{xo) / 0, then it 
is classical that / is strictly monotone. We can then write Bt = f~^{Xt) — rt. 
In particular, the random variables which are measurable with respect to Xt 
are measurable with respect to Bt and vice versa. On the other hand, by 
using a linear Gaussian regression, it is easy to show that D^Bt = Ht~^Bt 
(see also Proposition 5). Finally, the convergences (15) and (16) and the 
equality (14) allow us to conclude that we have (22). 
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Now, let us prove that the present is nondegenerate for X at t if and 
only if a{xo) 7^ 0. When (t(xo) = 0, it is clear that the present is degenerate 
at t (see the first part of this proof). On the other hand, if the present is 
degenerated at t, then there exists c G M such that 

Ht-^a o f{Bt + rt)Bt + rao f{Bt + rt) = c. 

By rearranging, we obtain that ao f{X){X -\-a) = (3 for some a, /3 G M, where 
X = Bt + rt. By using the fact that X has a strictly positive density on R, 
we deduce that a o f{x){x + a) = (3 for any x G M. Necessarily, /? = (with 
X = —a) and then f = a o f = 0. We deduce that / is constant and then 
that / = xo, that is, ct(xo) = 0. 

Finally, if Ht~^a{Xt)Bt + b{Xt) = a{Xt){Ht~^Bt + r) = a.s. for any i, 
then a{Xt) = = h{Xt) a.s. for any t and Xt = xq a.s. for any t; see (8). In 
other words, the collection of a-fields T = {Tf , t € (0,T)} is discriminat- 
ing. D 

Let us now describe a more general case. 

Theorem 20. Let X eE be given by (8) and let t G (0,T). Assume, 
moreover, that b £C^ and that a £C'^ is elliptic, that is, satisfies infa;gR |cr(x)| > 
0. Then X admits a forward and a backward stochastic derivative with re- 
spect to the present t, given by 

DlXt = D^Xt 



b{Xt) + H 



(T{Xt 



(23) 



Xt 



E 



dy 
o-(y) 

•*6 



Jo T 



ds 



t rt 



3H 



JO 

rt 



(3" (s) 5Br ds 



t I (3^{t)6Br 

Jo 



Xt 



where 



/?.(*)= Of/ 



b'a - ba' 



a 



{Xs)ls>.ds]it). 



Recall that Oh is defined by (A). 



Proof. First, note that j3^{t) belongs to the domain of the divergence 
operator 5^ , due to the additional hypothesis on b and a. We provide only 
the proof for D'^Xt, the computation for D_Xt being similar. Folmer [6], 
Section 4, tackles the problem of the computation of the time reversed drift 
of a non-Markovian diffusion by means of a Girsanov transformation and the 
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Malliavin calculus. Our proof uses such a strategy, coupled with the transfer 
principle. 

First step. Assume that a = 1. Using the transfer principle and the isom- 
etry K-h, it holds that 

rt 



where 



Here, we set 



Xt= f KH{t,s)dYs^ 
Jo 



Yt = Wt+ Car dr. 
Jo 



Kjl j'^b{Xs)dsyr). 



We know (see [13], Theorem 2) that the process X = (-^t)tGfo,Tl is a fractional 
Brownian motion under the new probability measure Q = G • P, where 

G = exp ( - / Os dWs — \j fls ds I . 

Using the integration by parts, of Malliavin calculus, we can write, for g : M ^> 
MgC^, 

E[(X,+;, - Xi)5(Xt)] = EQ[G-i5(^t)<5^(l[t,t+h])] 

= EQ[G-i(l[,,,+,],I)^5(^t))H] 

+ EQ[5(Xi)(l[i,i+^],Z)^G-i)^] 

= ^y {.^t)]{\t,t+h],\Q,t])H 

But K:%D^G~^ = D^G~^ (transfer principle). Since 



we have 



G = exp ( / as dYs — \ j a^ds 



G X dY{G~^) =at+ f Djas dV^ - f asD^as ds 
Jo Jo 

j-T 

= at+ DjasdWs. 
Jo 



Moreover 
rT 



[ DjardWr= [ {K.*HDfa){r)dWr= I L»f o^<5S,. := $(s) 
Jo Jo Jo 
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and 

(^ffl[o,t])(s) = ^i/(i,s)l[o,t](s). 
Therefore, 

{^*H'^[t,t+h],G)CHD' G~ )^2 

= {JCHa){t + h)- (ICnam + {JCH'^){t + h) - (/CH$)(t) 

= J biX^)du + {ICH<^){t + h)-ilCH^){t). 

By the stochastic Fubini theorem, we have {OH^){t) = /q {OHD.aj.){t) 53^- 
We set 

/?f (i) = {OnD^arm = (ohJ^ 6'(X,)l,>.(is)(t). 

We then deduce that 

E[iXt+f,-Xt)g{Xt)] 

(24) =E[g'{Xt)]m^t+^,l[o,t])H 

ft+h 



E 



9{Xt)U^ b{Xs)ds + J^ J P^{s)6Brds 



By developing E[Xt g{Xt)] as in (24), we obtain 
t2«E[5'(X,)]=E 

Then 

E[AhXt\Xt] 



g{Xt)[Xt- f b{Xs)ds- f f p^{s)6Brds 
\ Jo Jo Jo 



t rT 



h~\^^t+h]Am)n[^t-E ^ b{Xs)ds + J^ ^ P^{s)6Brds 



Xt 



t+h 



t+h rT 



+ h-^E I ' b{Xs)ds+ I / (5^{s)5BrdsXt 
Jt Jt Jo 

We deduce that E[/S.hXt \ Xt\ converges in probabihty, as /i J, 0, to 

t rt rT __ rT 



b{Xt) + jXt-jE 



b{Xs)ds+ I I /3^{s)6Brds- f P^{t)5B, 
Jo Jo Jo 



Xt 



Since hm/j|oE[A/iXf | Xt] does not depend on T, we finahy obtain (23) in 
the particular case where a = 1 by letting T J, i. 

Second step. Assume that a does not vanish. Set Yj = h{Xt), where 
h{x) = Jq -^T^- Using the change of variables formula, we obtain that Y 



satisfies 



Yt = yo + Bt+ [ -oh-\Ys)ds, te[0,T]. 

Jo CF 
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Since, on the one hand, the cj-fields generated by Xt and Yj are the same 
and, on the other hand, X has a-Holder continuous paths with a > 1/2, we 
have 

DlXt = a{Xt)DlYt. 

Expression (23) is then a consequence of the first step of the proof. D 

Remark 21. When a does not vanish and b = ra with r G M, we can 
apply either Proposition 19 or Theorem 20 to compute D^Xt. Of course, the 
conclusions are the same. Indeed, since we have, in this case, b'a — ha' = 
and /(f ' J^ = Bt + rt [since Xt = f{Bt + rt) with / satisfying /' = a o /], 

formula (23) can be simplified to (22). 

Compared to the case where a and h are proportional, here, it is more 
difficult to decide if the present (i.e., the collection of cr-fields generated by 
Xt) is discriminating or not. 

In the framework of the stochastic embedding of dynamical systems in- 
troduced in [2], the set of processes, called Nelson dijferentiable processes, 
which satisfy the equality between stochastic forward and stochastic back- 
ward derivatives plays a fundamental role (see [3], Chapters 3 and 7). We 
stress on the fact that solutions of stochastic differential equations driven 
by a fractional Brownian motion with H > 1/2 provide examples of Nelson 
differentiable processes which are not absolutely continuous. 
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